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ABSTRACT 

We describe the results of a new, wide-field survey for z = 3.1 Lya emission- 
line galaxies (LAEs) in the Extended Chandra Deep Field South (ECDF-S). By 
using a nearly top-hat 5010 A filter and complementary broadband photome- 
try from the MUSYC survey, we identify a complete sample of 141 objects with 
monochromatic fluxes brighter than 2.4 x 10~^^ ergs cm~^ and observers- frame 
equivalent widths greater than ~ 80 A (i.e., 20 A in the rest- frame of Lyct). The 
bright-end of this dataset is dominated by x-ray sources and foreground objects 
with GALEX detections, but when these interlopers are removed, we are still left 
with a sample of 130 LAE candidates, 39 of which have spectroscopic confirma- 
tions. This sample overlaps the set of objects found in an earlier ECDF-S survey, 
but due to our filter's redder bandpass, it also includes 68 previously uncataloged 
sources. We confirm earlier measurements of the z = 3.1 LAE emission-line lumi- 
nosity function, and show that an apparent anti-correlation between equivalent 
width and continuum brightness is likely due to the effect of correlated errors 
in our heteroskedastic dataset. Finally, we compare the properties oi z — 3.1 
LAEs to LAEs found at z — 2.1. We show that in the ~ 1 Gyr after 2; ~ 3, 
the LAE luminosity function evolved significantly, with L* fading by ~ 0.4 mag, 
the number density of sources with L > 1.5 x 10^^ ergs s~^ declining by ~ 50%, 
and the equivalent width scale-length contracting from 70^5 A to 50lg A. When 
combined with literature results, our observations demonstrate that over the red- 
shift range ^ ~ to 2; ~ 4, LAEs contain less than ~ 10% of the star-formation 
rate density of the universe. 

Subject headings: cosmology: observations; galaxies: formation; galaxies: high- 
redshift; galaxies: luminosity function, mass function 
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Introduction 



I n the local univers e, relatively few galaxies have Lya present in emission (jPeharveng et al. 



20081 : ICowie et al.ll2010l ). This is simply due to the nature of the line: since Lya involves the 
ground state of hydrogen, each photon must resonantly scatter many times before escaping 
into intergalactic space. As a result, even a small amount of dust can quench the emission 
leaving the galaxy. 

In the high-redshift universe, however, the presence of Lyct emission is common. Sur- 
veys of Lyman break galaxies (LBGs) have demon strated that roughly half of these bright , 
high- mass systems have Lya visible in emission ( IShapley et al.l 120031 : iReddy et al.l |2008| ) . 
Moreover, lower-mass Lya em itting objects with ex t remely high mass- specifi c star-forming 
rates exist in large numbers ( IGawiser et al.l l2006ai 120071 : lOuchi et al.l l2008l ). These Lya 
emitters (LAEs) are an important constituent of the high-redshift universe. Not only do 
LAEs sample a wider range of the galaxy luminosity function than LBGs, but they are 
also more amenable to spectroscopy. Consequently, LAEs are powerful probes of large scale 
structure, and perhaps the best tracers we have for the signature of baryonic acoustic oscilla- 
tions. In fact, LA Es are the targets of the Hobby-Eb erly Telescope Dark E nergy Experiment 
( iHill et al.ll2008l ). a Stage III Dark Energy project ( lAlbrecht et al.l 120061 ) which will survey 
the universe between 1.9 < z < 3.5. 

Unfortunately, we do not yet have a clear understanding of the global properties of 
LAEs. Although their space density, luminosity function, and equivalent width distribution 
are known at redshifts z > 3.1, the evolution of these properties towards lower redshifts is 
uncertain. For instance, the res ults of a intermediat e-bandpass (FWHM = 130 A) survey 
with the ESO 2.2-m telescope ( iNilsson et al.l l2009a( ) imply that LAEs become rarer and 
redder by z ~ 2.3, but the d ouble-blind VLT observations of Lya and Ha emiss ion argue 
otherwise ( iHayes et al.ll2010l ). Similarly, while the analysis of iNilsson et al.l (120111 ) suggests 
that the mean stellar mass of LAEs become s greater at lower redshift, a c omparison of 
narrow -band selected LAE samples at 2; = 3.1 (IGawiser et al.ll2007l ) and z = 2.l lGuaita et al. 
( I2OIII ) show no evolution in mass. 



Here we present the results of a new narrow-band survey for z = 3.1 LAEs in the 
Extended Chandra Deep Field South (ECDF-S). In §2, we describe our observations, review 
the techniques used to detect emission-line galaxies, and give a list of LAE candidates, 130 
of which are members of a statistically complete sample with Lya fluxes greater than 2.4 x 
10"^'' ergs cm~^ s~^ and rest-frame equivalent widths greater than 20 A. In §3, we present 
the z ^ 3.1 LAE luminosity function and show the distribution of Lya equivalent widths as 
a function of line luminosity and continuum lumin osity. In ^4, we comp are these functions to 
those of the sample of 2; ~ 2.1 objects obtained by lGuaita et al.l ( 20101 ). We show that there 
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has been significant evolution betwee n the two epochs in L*, the characteristic emission- 
hne luminosity of the ISchechterl ( 119761 ) function, 0, the integrated number density of bright 
LAEs, and Wq, the e- folding scale length of the rest-frame equivalent width distribution. In 
§5, we examine the LAEs' star-formation rates, using both their rest-frame UV emission and 
Lya, and in §6, we compare the star-formation rate density of the LAE population with 
that of the universe as a whole. We conclude by summarizing our results and arguing that 
LAEs represent the extreme low-metallicity, low-extinction, faint-end tail of the high-redshift 
star-forming galaxy population. 



For this paper, we assume a ACDM cosmology, with 
70 km s"^ Mpc-^ 



0.7, Qm = 0.3 and Hq 



2. The z = 3.1 Sample 



The luminosity function and equivalent w'idth distribution of z = 3.1 LAEs is fairly 
well-known f rom t he narrow-band surveys of iGronwall et al.l (120071 . hereafter Gr07) and 
Ouchi et al.l (120081 ). However these observations, which have identified more than 200 LAE 
candidates with little contamination from foreground objects, suffer from larger than normal 
photometric uncertainties due to the Gaussian-like transmission profiles of their narrow-band 
filters. To confirm these results, and to increase the survey sample at z ~ 3.1, we extended 
the Gr07 study by re-imaging the ECDF-S, this time with a 57 A full-width-half-maximum 
(FWHM) nearly top-hat filter centered near 5010 A. As with Gr07, the data were collected 
with the MOSAIC 11 CCD camera on the CTIO Blanco 4-m telescope, and consisted of a 
set of 47 images taken in I'.'l seeing for a total narrow-band exposure time of 15.67 hr. A 
log of the observations is given in Table [H a figure showing our filter's transmission curve, 
compared to that used by Gr07, appears in Figure [H 

The procedures used to reduce the data, identify line emitters, a nd measure LAE bright- 
nesses were identical to those used by Gr07 and discussed in detail by lFeldmeier et al.l (120031 ). 
After de-biasing, fiat-fielding, and aligning the data, we co-added our narrow-band frames to 
create a master image that was clipped of cosmic rays. This frame w as then compared to a 
deep B+V continuum image provided by the MUSYC collaboration (iGawiser et al.ll2006bl ). 
and the routines of DAOPHOT were used to create a color-magnitude diagram of all sources 
found on the narrow-band image. Objects with narrow-band minus continuum colors less 
than 0.93 mag in the AB system were flagged as possible emission-line sources (see Figure [2]). 
Then, to detect those LAEs whose photometry was compromised by object blending, sky 
gradients, and/or source confusion, we re-ran our detection algorithm on a "difference" frame 
made by subtracting a scaled version of the B+V continuum image from the narrow-band 
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frame. In this case, our detection algorithm was set to flag all objects brighter than four 
times the standard deviation of the local sky. This difference technique produced a sample 
of several thousand additional candidates, most of which were below our equivalent width 
threshold. However, when the weak emission-line sources were excluded, the result was a 
sample of objects that was ~ 15% larger than that produced by the two-color method alone. 

For both detection algorithms, we intentionally biased our parameters to identify faint 
sources at the expense of false detections (i.e., we introduced Type I errors into our dataset). 
To compensate for this, we visually inspected each candidate, both on the co-added narrow- 
band and broadband frames, and on sub-samples of the original images. This step excluded 
most of the false positives, which were typically at the frame limit, and artifacts produced 
by cross-talk associated with the CCD electronics. 

To derive the line fluxes and equivalent widths of our LAE candidates, we took advan- 
tage of the fact that Gr07 measured the 5000 A AB magnitudes of several ECDF-S field stars 
by comparing t heir large aperture instrumental magnitudes to those of spectrophotometric 



standard stars (IStond 119771 : iHamuy et al.l Il992l ) taken throughout their survey. To place 



our new 5010 A narrow-band observations on the same photometric system, we therefore 
measured our LAE brightnesses with respect to these same stars. We then converted our 
narrow-band 5010 AB magnitudes into monochromatic fluxes (-F5010) by comparing the fil- 
ter's integral transmission (which is the relevant quantity for field star photometry) to it s 



monochromatic transmission at the center of the bandpass (see Gr07 and lJacoby et al.lll987l ). 
Note that this simple procedure is only exact for top-hat filters where the monochromatic 
transmission is insensitive to wavelength: as Gr07 showed, observations through filters with 
a Gaussian-shaped transmission curve require a more sophisticated approach. 

Along with monochromatic flux, we also derived photometric equivalent widths for all 
our LAE candidates by comparing their narrow-band AB flux densities, /501O) to their B + V 
continuum flux densities, fs+v via 

ew=|^^-i1aa (1) 



I Jb+v 

where AA = 57 A approximates the contribution of the galaxy's underlying continuum within 
the narrow-band filter's bandpass. For consistency with Gr07, we then excluded all sources 
with equivalent widths less than 20 A in the rest-frame of Lya (approximately 80 A in the 
observers frame), and derived equatorial positions for the remaining obje cts with respect to 



reference stars in the USNO-Bl.O astrometric catalog (IMonet et al.l 120031 ). This left us with 
a sample of 199 LAE candidates. 

To measure the completeness of our sample as a function of magnitude, we followed 
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the procedures described by iFeldmeier et al.l (120031 ) and performed a series of artificial star 



experiments. By adding 1,000,000 stars (2000 at a time) to our narrow-band frame, and 
re-running our detection algorithms, we were able to compute the flux level below which the 
object recovery fraction dropped below 90% and 50%. The 90% value, which corresponds to 
an AB magnitude of msoio ~ 25.05, and a monochromatic flux of 2.4 x 10"^'' ergs cm~^ s^^ 
(logFsoio = —16.61) is our limiting flux for statistical completeness. A total of 141 objects 
are brighter than this limit. Those emission-line objects satisfying our 80 A equivalent width 
threshold and brighter than our completeness limit are listed in Table El additional objects 
below the completeness limit are given in Table [31 Our estimated photometric errors, as a 
function of log flux, appear in Table HI 

To confirm the robustness of our measurements, and to double-check the transmission 
curve of our narrow-band filter, we compared our sample of Lya emitters to the list of LAEs 
published by Gr07. As Figure [H shows, the volumes studied in these two surveys have a 
great deal of overlap: a quick inspection of the two curves suggests that ~ 40% of the 
LAEs found in our new survey should appear in the Gr07 tables. Indeed, this is the case: 
of the 199 candidates detected, 80 are present in the previous list of 2; ~ 3.1 emission-line 
sources. Moreover, 39 of the objects common to both surveys already have spectroscopic 
confirmations. In this restricted sample, the lowest redshift object has its emission line at 
4981 A (z = 3.096), where the transmission of our filter falls to roughly 25% of its peak 
value. Conversely, our highest redshift source has Lya at 5023 A (z = 3.131), i.e., where the 
Gr07 filter transmission is ~ 20% of its peak. A comparison of the inferred monochromatic 
fluxes for this subset of objects, where the wavelength of the Lya emission-line is known, 
also shows excellent agreement (see Figure [3l). This consistency argues that the transmission 
curves displayed in Figure [H are accurate, and the photometric errors tabulated in Table HI 
are reasonable. 



2.1. Chandra, GALEX, and Swift Detections 



Figure HI shows the flux distribution of all our high-equivalent width (EW > 80 A) ob- 
jects. The most obvious feature is the presence of a small number of sources with monochro- 
matic fluxes brighter than ~ 10~^^ ergs cm~^ s~\ T his luminous tail is not nearly as 
pronounced in the z = 3.1 survey of lOuchi et al.l ( l2008l ). and not visible at all in the Gr07 
dataset. A closer examination, however, reveals that many of these extremely bright sources 
are foreground interlopers. 



To show this, we first removed the AGN from our sample by cross- correlating the po- 
sitions of our LAE candidates with those of the x-ray point sources found in the 2 Msec 
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exposure of the Chandra Deep Field South (ILuo et al.ll2008l) and the four 250 ksec ex posures 



of the Extended Chandra Deep Field South (iLehmer et al.l l2005l : IVirani et al.l |2006[ ) . Eight 
spatial coincidences were found, including one for our second brightes t source (a. z = 1.609 



AGN which entered our sample due to its strong C III] A 1909 emission; iBalestra et al.ll2010l ) 
and our third brightest source (which, at z = 0.977, must have been detected d ue to a change 



in br ightness between the time of our narrow-band and continuum exposures; ISzokoly et al. 



2004J ). In fact, at flux levels brighter than logF > —15.9, half of our LAE candidates (and 
the candidates of Gr07) are co-incident with x-ray identified AGN. These contaminants are 
noted in Table [2l 

We next compared our list of LAE candidates to the c atalogs of UV sources d etected by 
the GALEX JSchiminovich et al.ll2003[ ) and Swift fUVOT: Eoversten et all 120091 ) satellites. 
Both instruments have conducted deep (m^s ~ 25) surveys of the ECDF-S using near-UV 
(and, for GALEX, far-UV) filters. Since at z ~ 3.1, the near-UV corresponds to rest-frame 
wavelengths far beyond the Lyman break, any LAE candidate listed in these catalogs is most 
likely a foreground contaminant. Our cross-correlation procedure (which used a matching 
radius of 2"), produced ten coincidences, including five associated with the x-ray sources 
described above. Of the remaining objects, the most prominent is our LAE candidate 7^1, 
a marginally-resolved galaxy which is nearly a magnitude brighter than any other object in 



our sa mple. Fortuitously, this object already has a redshift: spectroscopy by IVanzella et al. 



(120081 ) confirms the source to be a z = 0.337 galaxy with extremely strong [O II] A3727 
emission (rest-frame equivalent width EWq ~ 73 A). Although such objects are rare, and 
none were found in the GrO 7 survey, a few are present in the magnitude-limited set of objects 
defined by iHogg et al.l (119981 ). The other four UV matches include two likely superpositions 
(i.e., objects containing a more-likely counterpart within the error circle), and two sources 
with estimated observers- frame equiv alent widths (82 A) barely above our threshold. (The 
spectroscopy of iLe Fevre et al.l (120051 ) confirms that one of them is foreground.) We identify 
each of these UV-bright sources in Table [21 but in the analysis that follows, we exclude only 
the x-ray sources and the three GALEX/Swift-UVOT objects with unambiguously identified 
counterparts. This leaves us with a statistically complete sample of 130 objects, 39 of which 
have spectroscopic confirmations. 



2.2. X-ray Stacking 

By cross-correlating our LAE positions with the x-ray catalog from Chandra, we were 
able to eliminate bright AGN from our candidate list. This, of course, does not exclude the 
possibility that many of our LAEs have low-level nuclear activity. To place a constraint on 
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this activity, we per formed a stacking analysis on the x-ray data, similar to that described by 



Brandt et al.l (120011 ) . We combined our new LAE dataset with that found by Gr07, excluded 



the known x-ray sources, and stacked the x-ray signal of all 233 LAEs in the statistically 
complete samples, using the new 4 Msec exposure of the Chandra Deep Field South and the 
250 ksec exposures of the surrounding region. The mean LAE of the resulting ~ 480 Msec 
effective exposure has an x-ray flux less than 7.24 x 10~^ counts s~^ in the soft band and 
1.23 X 10"^ counts s~^ in the hard band (Icr upper limits), implying luminosity limits of 
Lx < 5.3 X 10^° ergs s~^ and Lx < 2.5 x 10^^ ergs s~^, respectively. Clearly, few, if any of 
our LAEs harbor low-luminosity AGN. We also used these data to place a limit on the mean 
star-formation rate of the LAEs in our sample. Since z ~ 3.1 LAEs all have extremely high 
mass-specific star-formation rates (iGawiser et al.ll2007| ). the contributions of low-mass x-ray 
binaries to the total 2-10 keV emission of these systems should be negligible. This being the 
case, our upper limit on the hard x-ray flux from the LAE populatioi i yields a star-format ion 
rate of < 33 Mf;^ yr ~^ or < 53 Mq yr~^, depending on whether we use lLehmer et al.l (120101 ) or 



Persic et al.l ( 120041 ) for the calibration. Although this value is still several times larger than 



the star- formation rates derived from the LAE's rest-frame UV and Lya emission (see §5), 
the numbers do provide hope that with additional observations, we will be able to detect 
the x-ray emission from z ~ 3.1 LAE starbursts. 



3. The z = LAE Luminosity Function and Equivalent Width Distribution 



Figure [5] displays the distribution of continuum-subtracted Lya emission-line luminosi- 
ties for our new sample of z = 3.1 LAEs. To find the best-fit ISchechten (119761 ) function 
to these data, we followed the procedures described in detail by Gr07. First, we computed 
two convolution kernels, one (a Gaussian) representing the photometric uncertainty of our 
measurements (as a function of magnitude), and the other determined using the derivative 
of the inverse of the filter transmission function. (This latter kernel does not normalize to 
one, since it also includes the effect of equivalent width censorship.) We then applied these 
two kernels to a series of functions of the form 



(t){L/L*)d{L/L*) = (P* {L/L*f e-^'^*d{L/V), 



(2) 



treated each curve as a probability distribution, and computed the likelihood that the ob- 
served sample of Lya luminosities could be drawn from that distribution. Figure [S] displays 
our best-fit function, both before (the dashed line) and after the equivalent-width censorship 
(solid red line) . Also plott e d are the best-fit curves for the z = 3.1 LAEs found by Gr07 
(in green) and lOuchi et al.l (l2008l . in blue). The parameters of these functions ar e listed in 
Table [51 Note that for the current study, we have fixed the faint-end slope of our ISchechter 
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( 119761 ) function to a = —1.65. This is the most-hkely value inferred from the Gr07 data, 
once continuum subtraction has been apphed to their fluxes, a r id it is also consistent with 
the recent 2 < z < 6.6 LAE measurements of ICassata et al.l (120111 ). Since our data are 
~ 0.5 mag shallower than that of Gr07, we adopt their estimate for the faint-end slope. 

From the figure, it is immediately obvious that all three luminosity functions are similar, 
both in shape and in normalization, and all fit th e data extrernely we ll. Our most-likely value 
of log L* = 42.76±0.10 is similar to that found by lOuchi et al.l ( l2008l ) and less than ~ 0.1 mag 
brighter than that derived Gr07. Similarly, all three functions imply ~ 7x 10~^ Mpc~^ for the 
density of galaxies with rest-frame equivalent widths greater than 20 A and Lya fluxes greater 
than 2.4 x 10"^'^ ergs cm~^ s~^. Since our new dataset has substantial overlap with t hat o f 



Gr07, the consistency between these two samples is not surprising. The lOuchi et al.l (|2008[ ) 
survey, however, is an independent measurement. Even after we follow their recommendation 
and increase 0* by ~ 10% to compensate for their higher equivalent width limit, there is 
still excellent agreement in the functions. This is indicative of the robustness of the results. 



3.1. Equivalent Width versus Luminosity 



There has been some discussion in the literature concerning the relationship between 
Lyg equiva l ent w idth and continuum luminosity in high-redshift galaxies. Starting with 
Ando et al.l ( 120061 ) . a number of authors have claimed a deficit of contin uum bright, high- 
equivalent width objects, both in samp les of Lyman-break galaxies (e.g . . lAndo et al.ll2007l: 



Iwata et al.l l2007l: IVanzella et al.ll2009l) and Lya emitters ( lOuchi et al. 



2008 



Shioya et al. 



20091 ). However, iNilsson et al.l ( l2009bl ). via a series of Monte Carlo simulations, have argued 
that the observed correlations are not statistically significant, and there is no need to invoke 



the presence of exotic stellar populatio ns or a clumpy interstellar medium ( lOuchi et al.ll2008 



Mao et al.ll2007l : lKobayashi et al. I boiol ). We illustrate the problem in Figure |6l The left panel 
of the figure shows the LAE equivalent width measurements as a function of (continuum 
subtracted) Lya luminosity, for both the present sample and the sample of objects found by 
Gr07. The figure displays no evidence for a shift in the eq uivalent width distributi on as a 
function of luminosity. A 2-sample Anderson-Darling test ( IScholz fc Stephenslll987l ) which 



compares the equivalent widths of bright (L > 3 x 10^^ ergs s ^) and faint (1.25 x 10^^ < 
L < 3 X 10^^ ergs s~^) LAEs confirms that there is no significant difference between the two 



samples. Moreover, if we apply the non-parametric rank-order lEfron fc PetrosianI (119921 ) 
test to the data, we find no reason to reject the hypothesis that equivalent width and Lya 
luminosity are statistically independent quantities. 



In contrast, the middle panel of Figure [HI displays the distribution of Lya equivalent 
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widths as a function of continuum magnitude. The difference between the two diagrams is 
striking. While the lack of low equivalent-width sources with faint continuum magnitudes 
is simply a selection effect (the Lya emission-line is too faint to be seen), the absence of 
high equivalent- width, bright co ntinuum objects is real. If we apply the non-parametric 



rank-order statistic described by lEfron fc PetrosianI (119921 ) to our truncated dataset, then 



we find that the null hypothesis of independence between equivalent width and continuum 
magnitude is excluded at the 5a level (for the Gr07 sample) and the 7a level (for our new 
z = 3.1 dataset). In other words, equivalent width shows no correlation when plotted against 
line brightness, but a clear anti-correlation when compared to continuum brightness. 

The explanation for these seemingly contradictory results probably lies in the fact that 
equivalent width is a derived quantity formed from two photometric measurements, one for 
the emission line, and one for the continuum. Consequently, the points displayed in the 
first two panels of Figure [6] possess correlated errors. In the left-hand diagram, emission-line 
strength and equivalent width are positively correlated, so photometric errors will preferen- 
tially scatter objects into the high-luminosity, high-equivalent width region of the diagram. 
Conversely, in the middle panel, continuum strength and equivalent width are negatively 
correlated, so errors in the abscissa will move objects away from this portion of the figure. 
When this dichotomy is coupled with the heteroskedastic nature of the dataset - the largest 
uncertainties are associated with the objects having the faintest continuum magnitudes and 
therefore the highest equivalent widths - the result can be an apparent correlation in one 
figure, and a lack of correlation in another. 

The best way to investigate the systematics of equivalent width is to plot the two 
independent quantities in the relation - line strength and continuum strength - against one 
another. This is done in the right-hand panel of Figure [61 Here, there is no evidence of 
any systematic behavior with luminosity: the LAEs occupy a strip in the diagram whose 
thickness appears constant with luminosity. This argues strongly against the presence of a 
systematic dependence of equivalent width on luminosity. 



4. Evolution from 2; ~ 3 to 2; ~ 2 

To investigate the evolution of Lya emitters in the redshift range 2.1 < z < 3.1, we 
can compare the pr operties of our LAE candidates to those of the z ^ 2.1 objects found by 



Guaita et al.l (120101 . hereafter GulO) in their deep narrow-band survey of the ECDF-S. Like 
the observations described here, GulO used a 50 A wide FWHM filter (Ac = 3727) and the 
prime focus Mosaic II CCD camera of the CTIO 4-m telescope to identify a complete sample 
of candidate LAEs over ~ 1000 arcmin^ region of the sky. Also, like our 2; ~ 3.1 sample, the 
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GulO objects are defi ned via their color wi th respect to deep continuum images provided by 
the MUSYC survey (iGawiser et al.ll2006bl ). have rest-frame Lya equivalent widths greater 
than 20 A, and have minimum Lya luminosities of ~ 1.3 x 10^^ ergs s~^. Thus, in area, 
number, and depth, the GulO dataset is well-suited for a comparative study. 

GulO have already used the Chandra and GALEX source lists to eliminate some high- 
equivalent width interlopers from their sample: as with our 2; ~ 3.1 candidates, true z ~ 2.1 
galaxies without AGN should not only be undetected in the x-ray, but in the ultraviolet as 
well, as nearly all their near-UV fiux originates at wavelengths beyond the Lyman break. 
To these criteria, we added two others. First, we eliminated one LAE c andidate that wa s 
previously identified as a z = 2.0 QSO in the COMBO-17 source catalog JWolf et al.ll2004l ). 
EIS J033330.60-274819.3 is located at the extreme edge of the ECDF-S Field 1 x-ray point- 
ing, ~ 9' from the field center. Although the object is not cataloged as an x-ray source, it 
is (barely) detected in the hard x-ray band, and its im plied (rest-frarae) 2 keV to 2500 A 
fiux ratio is consistent with that expected for an AGN (IStrateva et al.l 120051 ). (Using a 1'.'5 
aperture, we derive Oox = "l-S^oi-) Moreover, the candidate in question is exceptionally 
luminous for an LAE, as it has a continuum-subtracted narrow-band fiux that is over a 
magnitude brighter than the next most luminous source in the GulO object list. Thus, even 
if the candidate is not an AGN, its extreme luminosity places it in a differ ent category frora 
all other candidates in the GulO sample, and unfittable via a traditional ISchechterl (Il976[ ) 
function. 

In addition to excluding this one object, we also examin ed each candidates' ap pear- 
ance on HST V^-band images taken as part of the GOODS (iGiavalisco et al.l 120041 ) and 
GEMS (IRix et al.l 120041 ) surveys. Since LAEs are, in general, much fainter in the continuum 
than Lyman-break galaxies, we might expect the sample of 2 ~ 2.1 objects to be smaller 
(in physical extent) than the l argest Lyman-break galaxies found by iFerguson et al.l (l2004r) 
and 



Ravindranath et al.l ( l2006l ). In fact, the HST-based morphological study of 



Bond et a 



( I2OIII ) classified these "large" objects as photometric interlopers, and iGuaita et al.l (120 111 ) 
used this discriminant in their analysis of the LAEs' spectral energy distributions. Thirty- 
four of the GulO objects fall into this category, but only four are among the brightest 100 
objects in the sample, and only one is more luminous than the sample's 90% completeness 
limit. Consequently, their inclusion or exclusion makes very little difference to the analysis. 
Below, we will present results both including and excluding these objects. 

Before fitting the luminosity function, we made one other modification to the GulO 
dataset. For consistency with our 2 ~ 3.1 photometry and with the analysis of Gr07, we 
chose to adopt GulO's point-source aperture fiuxes in our analysis, rather t han their cor- 
rected aperture fiuxes, which were derived using the photometric techniques of iGawiser et al. 
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( l2006bl ). While this latter procedure is most useful for selecting LAEs out of the general 
galaxy population, the former method is better suited for measuring the brightnesses of ex- 
ceedingly compact (< 0'.'25) sources on groundbased images. (Since the seeing on the GulO 
frames is roug hly twice the size of the largest LAE counterpart identified on HST frames 
of the region (IBond et al.l 120111 ). this change is certainly justified.) This technical revision 
reduces the derived magnitudes of a significant fraction of the GulO LAEs, sometimes by 
as much as ~ 0.3 mag. Our final 2; ~ 2.1 sample then consists of 37 LAEs brighter than 
the GulO 90% completeness limit of rriAB = 24.8, and 73 sources belonging to a statistically 
complete sample of objects with monochromatic, continuum-subtracted fluxes greater than 
2 X 10"-'^'' ergs cm~^ s"-*^. 

In Figu re [71 we fi t the GulO LAE monochromatic fluxes (with the "large" objects ex- 
cluded) to a lSchechterl (119761 ) function in exactly the same manner as for the 2; ~ 3.1 samples, 
setting a = —1.65, and using convolution kernels appropriate for the GulO photometric er- 
rors and filter transmission function. The first feature to note is the offset between the dashed 
line, which displays the best fit to the data, and the solid black line, which represents the 
implied luminosity function. This difference, which is almost 3 times larger than in Figure \5\ 
is due entirely to the properties of the narrow-band filter employed in the survey. Like the 
Gr07 observations, the GulO study was performed through a filter whose transmission curve 
is more Gaussian-shaped than top-hat. Consequently, LAEs with redshifts that place their 
emission-line off the center of the filter have their monochromatic fluxes (and equivalent 
widths) systematically underestimated. The result is an effective survey volume that is less 
than that implied by FWHM of the filter. (See Gr07 for a full description of the effect.) 

Even after we correct for this offset, it is apparent that there are fewer bright LAEs 
at 2; ~ 2.1 than predicted by the z ~ 3.1 luminosity functions. Our z ~ 2.1 value of 
logL* = 42.33 ± 0.12, is ~ 0.4 mag fainter than that measured at the higher redshift, and 
there are ~ 50% fewer LAEs with emission-line luminosities brighter than 1.5 x 10^^ ergs s~^. 
Even if the "large" LAE candidates are included in the fit, the result does not change: the in- 
clusion of these objects brightens L* by only 0.01 mag, and increases the normalization of the 
function by only 4%. The dramatic change in the luminosity function between z = 3.1 and 
z = 2.1 is also illustrated via th e likelihood c ontou rs of Figure El and is similar, though a bit 
larger, than the fall-off seen by iBlanc et al.l (1201 if ) in their blind integral-field spectroscopic 
survey of the COSM OS, GOODS - N, M UNICS-S2, and XMM-LSS fields. We note, how- 
ever, that though the iBlanc et al.l (120111 ) observations with the McDonald 2.7-m telescope 
dwarf all other LAE surveys in terms of survey volume ~ 10^ Mpc^), their measurement of 
LAE evolution is limited by the relatively small number of objects detected (~ 100 LAEs 
distributed roughly evenly between 1.9 < z < 3.8), and a very bright limiting luminosity 
(~ 10^^ ergs s~^). Similarly, our results are consistent with the z = 2.25 narrow-band pho- 
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tometry of iNilsson et al.l ( l2009al ). though again, the size of the telescope (2.2 m) and width 
of the interference filter (130 A FWHM) restricted their study to relatively luminous objects. 
Still, this survey did detect ~ 190 LAE candidates with little expected contamination, and 
the global properties of the sample sho uld be robus t . The only apparent contradiction with 



our data involes the z = 2.2 analysis of iHayes et al.l (120101 ). who used a combination of Lya 



Ha, and deep broadband imaging (the latter for photo-z estimates) to derive a much slower 
rate of evolution at the bright end of the LAE luminosity function. However, this result is 
based on only 21 LAEs, none brighter than 2 x 10^^ ergs s^^. Thus, it is still possible that 
all the results are consistent. 

More evidence of LAE evolution comes from the distribution of Lya equivalent widths. 
Figure displays a histogram of rest-frame equivalent widths for the present z = 3.1 sample, 
the Gr07 z = 3.1 sample, and the GulO z = 2.1 dataset. The two z = 3.1 distributions are 
statistically identical: K-S and Anderson-Darling tests can discern no difference between the 
datasets, and an exponential fit to the distributions, using the same maximum-likelihood 
procedures as Gr07, yields a rest-frame scale-length of wq = 64^7° A for the LAEs of our 
new survey. This is consistent with the value of Wq = TGlg^ A found by Gr07, and, when the 
two datasets are combined, the resulting equivalent width scale length for z = 3.1 LAEs is 
70^5 A. This contrasts sharply with the best-fit value of Wq = 50lg A for the z ~ 2.1 LAE 
sample. Non-parametric statistics, such as those associated with the K-S and Anderson- 
Darling tests, prove that the samples differ with greater than 99% confidence, confirming 
that some amount of evolution has occurred in the Gyr between z = 3.1 and z = 2.1. 



5. Lya and UV Star Formation Rates 



The strength of the Lya emission line and brightness of the UV continuum are both 
functions of star formation. Specificall y, the as s umpti ons of an optically thick interstellar 
medium. Case B recombination, and a ISalpeterl (119551 ) universal initial mass function lead 
to a relationship between Lya luminosity (in ergs s~^) and star- formation rate. 



SFR(Lya) = 9.1 x 10~^^L(Lya) Mq yr"^ 



(3) 



( Kennicutt 


1998; 


Brocklehurst 


1971) 



stars undergoing a constant rate of star formation predict 

SFR(UV) = 1.4 X 10-28 Mq yr^^ 



(4) 



where (in erg s s ^ Hz ^) repre sents the population's mid-UV (1500 A to 2800 A) contin- 



uum luminosity (IKennicutt 



19981 ). Thus, in the absence of dust, galaxies which are optically 
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thick to Lyman radiation should have Lya luminosities that are well-correlated with their 
continuum flux. 

Figure HU] compares the two quantities. From the figure, it is easy to see that at both 
redshifts, there is a wide range of star-formation rate ratios, and that the scatter in the 
diagram is much larger than that produced by simple photometric errors. Objects above the 
one-to-one line are ea sily explained by the loss of Lyg photons, most likely to encounters 
with dust grains (e.g.. ICharlot fc Falllll993HHaves et al.ll2010l : lAtek et al.ll2009l ). In contrast, 
objects below the line may be systems for which the "continuous star formation approxi- 
mation" in population synthesis models is inappropriate (i.e., objects with star- formation 
timescales less than ~ 10*^ yr) or ones where the dust is preferentia l ly embedded within 
large neutral gas cloud s with steep density profiles (e.g., 
Finkelstein et al.lboogl ). 



Neufeld 1991; Hansen fc Oh 2006 



Interestingly, the distribution of star-formation rate ratios displayed in Figure [TO] shows 
no evidence for evolution between z ~ 2.1 and z ~ 3.1. Our new z = 3.1 survey is complete 
to a continuum-subtracted monochromatic luminosity limit of 1.2 x 10^^ ergs s~^. For LAEs 
brighter than this, the median ratio of the UV star-formation rate to the Lya star-formation 
rate is 1.09lJ^;^^, where the uncertainty is determined via a bootstrap analysis. For z ^ 2.1 
systems, this ratio is l.lll^o ^5. Moreover, as Figure [TT] illustrates, there is no difference in the 
shape of the underlying distribution: while a number of physical mechanisms might cause a 
systematic shift in the ratio of the two star formation rate indicators, an Anderson-Darling 
test can find no evidence for any evolution, even at the 2 a level. 

Nevertheless, the change in the equivalent width distribution, does point to a decrease 
in the efficiency of Lya emission. As described above, Lya line fiux and mid-UV continuum 
fiux are both driven by a galaxy's star-formation rate. Moreover, in the wavelength range 
1250 A < A < 2600 A, a star-forming galaxy's spectral e nergy distribution (SEP) is well-fit 



by a simple power-law. Fa oc A^, with slope /3 ~ —2.3 (ICalzetti et al.lll994 iMeurer et al. 



I999I ). S ince is rather insensitive to a starburst 's duration, metallicity, or initial mass 
function (ILeitherer et al.lll999l : iBouwens et al.ll2009l ). any systematic shift in a population's 
Lya to 1250 A continuum fiux ratio, i.e., in its Lya equivalent width, must be due to some 
other parameter. Geometry is one such variable, since an increase in the opening angle for 
Lya emission would cause a typical line-of-sight to receive fewer photons. However, purely 
geometrical effects cannot explain the fading of the LAE luminosity function, nor a decline 
in the star formation rate density of the LAE population (see next section). The most likely 
explanation for these changes is dust, which can decrease /3, and also reduce the efficiency 
of line em ission by dest r oying Lya photons. This interpretation is supported by the SED 
analysis of lGuaita et al.l ( 1201 ll ). who found that z = 2.1 LAEs have slightly greater reddening 
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than their z = 3.1 counterparts (iGawiser et aLl 120071 ). 



6. The Star Formation Rate Density of LAEs 



Narrow-ba nd imaging (Ouchi et al.ll2008l . l2010l . Gr0 7, GulO), along wi th shtless ( jPeharveng et ah 



20081 ). long-sht (ICassata et aLll201ll ). and integral- field ( iBlanc et al.ll201ll ) spectroscopy now 
yield a coherent picture of the evolution of the LAE population and their star-formation 
rate density. To illustrate this, we computed the total observed Lya luminosities of the 
z = 2.1 and z = 3.1 LAE populations by integrating their best-fit ISchechterl (119761 ) func- 
tions of Table down to zero. We then used equation to convert these luminosities into 
star-formation rate densities, and compared our results to those of other LAE measurements 
and from estimates made using the rest-frame UV continuum. This comparison is shown 
in Figure [121 As the figure indicates, at 2; > 5, the star-formation rate densities one ob- 
tains from Lya are identical to those derived from the rest-frame UV; at these redshifts, the 
assumption of an optically thick Case B interstellar medium, little extinction, and a Lya 
escape fraction near unity appears reasonable. However, at a lookback time of ~ 12 Gyr, 
there is a rapid transition, as at all subsequent epochs the ratio of LAE to rest-frame UV 
star- formation rate density is < 0.1. The decline in the Lya-based measurements between 
z ^ 3.1 and z ~ 2.1 (and, indeed, to 2; ~ 0) roughly follows that for the universe as a whole. 

Figure [TO] shows that, in a differential sense, dust has very little effect on the 2 < z < 4 
LAE population: on average, the star-formation rates inferred from Lya are similar to 
those derived from the rest-frame UV. The dispersion in the diagram is large, suggesting 
the presence of some reddening (and perhaps implying that many of the LAEs have star- 
formation timescales less than ~ 10^ yr), but overall, the Lya escape fraction from these 
objects appears to be close to unity. On the other hand. Figure [12] proves that globally, 
Lya-based measurements of star formation yield values that are more than an order of 
magnitude less than those found from the rest-frame UV. Clearly, the bulk of star formation 
is occuring outside of the LAE population. As confirmed by their colors and nebular emission, 
LAEs repre sent the extreme low-extinction, low-rn e tallicity tail of the sta r-forming galaxy 
population (IGawiser et al.ll2007l : iGuaita et al.ll201ll : iFinkelstein et al.ll201ll ). 



7. Conclusion 



Our 5010 A narrow-band survey of the Extended Chandra Deep Field South has pro- 
duced 199 high equivalent width (EW > 80 A) objects, of which 141 are members of a statis- 
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tically complete sample with monochromatic fluxes brighter than 2.4 x 10^^^ ergs cm~^ s~^. 
X-ray and UV observations with the Chandra, GALEX, and Swift satellites demonstrate 
that most of our brightest candidates are AGN or foreground [O 11] emitters, but when we 
exclude these sources, we are still left with a complete sample of 130 candidates that are 
likely to be Lya emitting galaxies at z ~ 3.12. This sample overlaps that deflned by Gr07 
usmg a 4990 A narrow-band fllter, but includes an additional 68 candidates. This brings the 
total number of likely z ^ 3.1 LAEs known in the ECDF-S to 360, with 223 being members 
of statistically complete samples. The sample now extends in redshift space from z = 3.06 
to z ^ 3.14 and includes 72 objects with spectroscopic conflrmations. 

Using these data, we have explored the luminosity function and equivalent width dis- 
tribution of z ~ 3.1 LAEs. Both functions are in excellent agreement with previous deter- 
minations: our value for the Schechter cutoff (logL* = 42.76 ± 0.10) and equivalent width 
e-folding scale (w n = 64^7° A) are statistically identical to values in the literature (Gr07, 
Ouchi et al.ll2008l ) and the normalization of the functions are also consistent to ~ 25%. We 
show that an apparent apparent deflcit of continuum-bright, high equivalent-width objects 
is most likely an artifact caused by correlated errors in our heteroskedastic dataset. When 
the independent parameters of line-luminosity and continuum luminosity are compared, the 
apparent correlation disappears. 

We have also examined the evolution of the LAE luminosity function towards lower 
redshifts by comparing our sample of z ~ 3.1 Lya emitters to the z ~ 2.1 LAEs found by 
GulO. We show that in the ~ 1 Gyr between the two epochs, the number of LAEs with 
line luminosities brighter than 1.5 x lO''^ ergs s~^ declined by ~ 50%, and logL* faded by 
~ 0.4 mag. Thus, it appea rs that the large d e crease in th e frequency and str ength of galactic 
Lya emission observed by iDeharveng et al.l fcOOSl ) and ICowie et al.l (120101 ) at z ~ 0.3 has 
begun by z ~ 2. Similarly, there are significantly fewer z ~ 2.1 LAEs with extremely high 
equivalent widths: between z ^ 3.1 and z ^ 2.1, the scale length of the LAE equivalent width 
distribution decreased from 70^5 A to 50t^\^ A. This decrease suggests that z ^ 2.1 LAEs 
have slightly more dust than their z ~ 3.1 cou nterparts. This hypo t hesis is supported b y the 
results of spectral energy distribution fitting ( iGawiser et al.l 120071 : iGuaita et al.l 120111 ) . but 
is not confirmed by an analysis of the objects' UV to Lya-based star formation rate ratios. 

Finally, by integrating our LAE luminosity functions and then combining our measure- 
ments with other data from the literature, we have estimated how the contribution of LAEs 
to the star-formation rate density of the universe changes with redshift. At z > 5, LAEs are 
responsible for most of the star-formation rate density, b ut in the redshift ra nge 2 < 2 < 3, 
this fraction has decreased to ~ 0.1. As pointed out by iBlanc et al.l (120 111 ), this suggests 
that globally, there has been strong evolution in the Lya escape fraction. This is consistent 
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with other results based on LAE photometry and spectroscopy: at the epochs considered 
here, Lya emitters represent the extreme low-metalhc i ty, low-extinction, f a int-end tail of the 



star-f orming galaxy population ( IGawiser et al.l 120071 : iGuaita et al.l 1201 ll : iFinkelstein et al. 

201 ih . 
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Fig. 1. — Transmission curve of our narrow-band A5010 filter i n the f/2.7 beani of the CTIO 
4-m telescope, compared to that of the 4990 A filter used by iGronwall et al.l (120071 ) . Our 
filter's nearly top-hat profile simplifies the analysis of the current survey. 
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Fig. 2. — Excess emission in the narrow-band A5010 filter over the continuum for objects 
in our survey field. The abscissa gives A5010 AB magnitude, while the ordinate shows the 
difference between this magnitude and the magnitude on a B + V continuum frame. Our 
narrow-band completeness limit of 2.4 x 10~^^ ergs cm~^ s~^ is represented by a vertical 
line; our equivalent width limit of 80 A is shown via the horizontal line at msoio — f^B+v = 
—0.93. The magenta curve shows the expected 3 a errors in the photometry. Candidate 
emission line galaxies are denoted as red circles; the green circle shows an x-ray source 
detected by Chandra, the dark blue circles denote ultraviolet sources cataloged by GALEX 
and/or Swift/UVOT, and the light blue circles show objects detected in both the x-ray (with 
Chandra) and UV (with GALEX and/or Swift/UVOT). 
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Fig. 3. — A comparison of the continuum-subtracted emission-line fluxes of 39 spectroscop- 
ically confirmed LAEs detected with both a narrow-band A4990 filter and our narrow-band 
A5010 filter. The scatter is consistent with the photometric precision of the measurements. 
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Fig. 4. — The monochromatic flux distribution of objects with observed equivalent widths 
EW > 80 A. The brightest source is a resolved [O II] emitter; the next three brightest objects 
are x-ray sources. The open circles represent data beyond our completeness limit. 
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Fig. 5. — The differential luminosity function of z = 3.1 Lya galaxies with rest frame 
equivalent widths greater than 20 A binned into 0.2 mag intervals. The points give the 
density of objects under the assumption that our filter's FWHM defines the survey volume; 
the open circles r epresent data be yond our completeness limit. The solid red curve shows 
our input best-fit ISchechten (Il976l ) luminosity function, while the dashed line illustrates the 
shape and normalization of this function after correcting for the effects of photometric error 
and censo ring by our filter's non-square transmission curve. The luminosity functions of 
Gr07 and lOuchi et al.l ( 120081 ) are also plotted, though the curve for the latter is based on a 
rest-frame equivalent width limit of 64 A. For consistency with our selection criteria, this 
curve should be increased by > 10%. 
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Fig. 6. — The left-hand panel displays the distribution of Lya rest-frame equivalent widths as 
a function of (continuum subtracted) Lya emission-line luminosity, the middle panel shows 
these same equivalent widths as a function of continuum magnitude, and the right-hand panel 
plots Lya luminosity versus continuum luminosity. Crosses are values from Gr07, squares 
are from the current study, and circles represent LAE candidates found in both surveys. The 
red points show spectroscopically confirmed objects. The lines show the completeness limits 
of the surveys (solid for this survey, dashed for Gr07), and the arrows in the left and middle 
panels show how the errors in equivalent width and luminosity are expected to correlate. The 
heteroskedastic nature of the uncertainties makes a simple interpretation of the equivalent 
width diagrams difficult, but the last panel strongly suggests that there is no correlation 
between any of the quantities. 
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Fig. 7. — The differential luminosity function of 2; = 2.1 L AEs with rest - frame equivalent 
widths greater than 20 A, inferred from the observations of iGuaita et al.l f 2010l ). The data 
have been binned into 0.2 mag intervals; points beyond the completeness limit, or where 
photometric errors may cause an excess amount of contamination are shown as open circles. 
The open circle at the bright end represents the one extremely luminous object tha t has been 
classifi ed as a 2; = 2 QSO by I Wolf et al.l (120041 1. The solid red line is the best-fit ISchechter 
(119761 ) function, while the dashed line shows the shape and normalization of this function 
after correcting for the effects of photometric error and censoring by our filter's non-square 
transmission curve. The shaded area represents the range of z ~ 3.1 luminosi t y f mi ctions, 
and the blue line is the Schechter fit to 2; ~ 2.2 LAEs found by iHayes et al.l (120101 ). The 
data suggest that L* has faded by ~ 0.4 mag between z ~ 3.1 and z ~ 2.1. 



- 31 - 



.002 



.0015 — 



.001 — 



A .0005 
9- 







1 — I — r 



1 — I — I — r 



1 — \ — r 




J I L 



J I I L 



J \ I L 



43.5 



43 



42.5 



42 



Log L (ergs s ^) 



Fig. 8. — A comparison of the likelihood confidence contours of ISchechten (119761 ) function 
fits to the observed distribution of Lja luminosities. The red contours, drawn at 0.5 a 
intervals, show the results for our z = 3.1 LAE sample; the blue curves display the same 
information for the GulO z = 2.1 dataset. The abscissa r epresents t h e con tinuum-subtracted 
value for logL*; the ordinate displays the integral of the ISchechterl ( 119761 ) function down to 



a luminosity of 1.5 x 10 ergs s . Between the two redshifts, the density of LAEs has 
decreased by more than 50%. 
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Rest Frame Equivalent Width (A) 

Fig. 9. — Equivalent width histograms for our new sample of z ~ 3.1 LAEs, compared to 
that measured by Gr07 (for 2; ~ 3.1) and GulO (for 2; ~ 2.1). The curves are the best-fitting 
exponentials, corrected for the effects of photometric error and censoring by the filters' non- 
square transmission curves. There are significantly fewer high equivalent width objects at 
z ~ 2.1 than at 2; ~ 3.1. 
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Fig. 10. — A comparison of the star formation rates derived from Lya emission (under 
Case B recombination) and the UV continuum at 1570 A. The left-hand panel shows z ~ 
2.1 LAEs from GulO, while the right-hand panel displays z ~ 3.1 objects, with crosses 
representing LAEs from Gr07, squares denoting the sources found in the current study, and 
circles showing candidates present in both surveys. The red points show LAEs that have been 
spectroscopically confirmed. The diagonal dashed line shows where the two measurements 
are equal. The flux limits are shown via the vertical blue line (in the case of the z = 3.1 
data, solid for this survey, dashed for Gr07); the approximate location of our equivalent 
width threshold is shown via the diagonal blue line. Note that, although the two indicators 
are correlated, there is a substantial amount of scatter in the relation. 
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Fig. 11. — The ratio of the star-formation rate derived from the UV continuum to that 
derived from Lja (under Case B recomination) for all LAEs with line luminosities brighter 
than logL = 42.08. The dark bins represent objects for which we only have upper limits; 
the dashed lines show the median values of the distributions: 1.09 for z = 3.1 LAEs and 1.11 
for the z = 2.1 objects. An Anderson-Darling test cannot detect any difference between the 
two distributions. 
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Fig. 12. — Measurements of the observed LAE star-formation rate density as a function 
of redshift. Each point is derived from by integrating the Schechter function down to zero 
luminosity, and for consistency, each assumes a f aint-end slope of a = —1.65. The shaded 
area displays the results from the rest-frame UV (iBouwens et al.ll2010l ). At z > 6, the LAE 
and UV-based measurements of the star-formation rate density are equal, but by z ~ 4, this 
ratio has transitioned to a value of < 0.1. 
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Table 1. Log of Narrow-band Observations 





Exposure 


Seeing 




UT Date 


(min) 


(") 


Notes 


15 Oct 2006 


9 X 20 


I'/l 


Through clouds 


17 Oct 2006 


8 X 20 


I'.'l 


Photometric 


19 Oct 2006 


11 X 20 


I'.'l 


Photometric 


21 Oct 2006 


12 X 20 


I'/l 


Photometric 


20 Nov 2006 


3 X 20 


I'/l 


Photometric 


22 Nov 2006 


4 X 20 


l'/2 


Light cirrus 
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Table 2. Candidate Lya Emitters: The Statistically Complete Sample 



ID 


q(2000) 


<5(2000) 


Log -F5010 


Observed EW (A) 


z 


Comments^ 


1 


3:32:17.57 


-27:49:40.9 


-15.227 


98 


0.337 


GALEX, UVOT source; Resolved 


2 


3:32:04.06 


-27:37:25.4 


-15.523 


124 


0.977 


x-ray, GALEX source 


3 


3:31:46.18 


-27:57:08.6 


-15.523 


129 




x-ray, GALEX, UVOT source 


4 


3:33:07.61 


-27:51:26.9 


-15.707 


153 


1.609 


x-ray, UVOT source; 


5 


3:31:40.14 


-28:03:07.5 


-15.723 


249 




G-009 


6 


3:31:34.73 


-27:56:21.8 


-15.743 


377 


3.131 


G-019 


7 


3:33:12.72 


-27:42:46.9 


-15.767 


389 




G-002 


8 


3:31:51.04 


-27:57:15.7 


-15.771 


195 




x-ray, GALEX source 


9 


3:33:04.86 


-27:45:54.9 


-15.811 


97 






10 


3:32:18.80 


-27:42:48.2 


-15.851 


230 


3.116 


G-004 



^G-XXX denotes identification number in lGronwall et aL I liool) 
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Table 3. Candidate Lya Emitters: Objects Below the Completeness Limit 



ID 


a(2000) 


5(2000) 


Log F5010 


Observed EW (A) 


z 


Comments^ 


142 


3:32:06.86 


-27:53:33.4 


-16.619 


611 






143 


3:32:58.57 


-27:41:31.0 


-16.619 


307 




G-221 


144 


3:31:52.83 


-27:45:18.6 


-16.623 


186 


3.106 


G-050 


145 


3:33:12.89 


-28:03:08.2 


-16.627 


739 






146 


3:33:18.56 


-27:33:49.5 


-16.635 


926 






147 


3:33:21.69 


-27:45:06.9 


-16.635 


163 






148 


3:33:02.99 


-28:02:34.9 


-16.635 


136 






149 


3:32:57.35 


-27:51:42.1 


-16.635 


331 




G-102 


150 


3:33:38.74 


-27:47:02.5 


-16.643 


460 




G-142 


151 


3:32:14.57 


-27:45:52.4 


-16.647 


314 


3.118 


G-090 



''G-XXX denotes identification number in lGronwall et aL 1 ||2007| ') 



- 39 - 



Table 4. Photometric Uncertainties 



Log F5010 


cr (mag) 


Log F5010 


a (mag) 


-15.50 


0.012 


-16.30 


0.079 


-15.60 


0.016 


-16.40 


0.095 


-15.70 


0.022 


-16.50 


0.115 


-15.80 


0.029 


-16.60 


0.134 


-15.90 


0.036 


-16.70 


0.167 


-16.00 


0.044 


-16.80 


0.198 


-16.10 


0.052 


-16.90 


0.235 


-16.20 


0.066 


-17.00 


0.288 



Table 5. Best-Fit Schechter Function Parameters 





Sample 


log L* 
(ergs s~^) 


a 


log0* 
(Mpc-3) 


(ergs s~^ Mpc~^) 


z 


= 3.095 (Gr07) 


42.70 ±0.10 


-1.65 


-3.14 ±0.04 


39.97 ±0.11 


z 


= 3.113 


42.76 ±0.10 


-1.65 


-3.17 ±0.05 


39.99 ±0.11 


z 


= 3.1 rOuchi et al. 2008)'^ 


42.82 ±0.06 


-1.65 


-3.15±0.11 


40.08 ±0.13 


z 


= 2.063 (GulO) 


42.33 ±0.12 


-1.65 


-2.86 ±0.05 


39.59 ±0.13 



Values are interpolated to a = —1.65, and appropriate for a 64 A rest-frame equiv- 
alent width limit. Following their suggestion, these numbers should be increased by 
~ 10% to approximate our 20 A limit. 



